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ABSTRACT 

We analyze the absolute magnitude (M^) and color [u — r) of low redshift 
{z < 0.06) galaxies in the Sloan Digital Sky Survey Data Release 6. Galaxies 
with nearly exponential profiles (Sloan parameter fracDeV < 0.1) all fall on the 
blue sequence of the color - magnitude diagram; if, in addition, these exponential 
galaxies have Mr < —19, they show a dependence of m — r color on apparent 
axis ratio q expected for a dusty disk galaxy. By fitting luminosity functions 
for exponential galaxies with different values of q, we find that the dimming 
is well described by the relation AM^ = l-27(logg)^, rather than the AM oc 
log q law that is frequently assumed. When the absolute magnitudes of bright 
exponential galaxies are corrected to their "face-on" value, M/ = — AM^, 
the average u — r color is linearly dependent on M/ for a given value of q. 
Nearly face-on exponential galaxies (g > 0.9) have a shallow dependence of mean 
u — r color on M/ (0.096 magnitudes redder for every magnitude brighter); by 
comparison, nearly edge-on exponential galaxies {q < 0.3) are 0.265 magnitudes 
redder for every magnitude brighter. When the dimming law AM^ oc (logg)^ 
is used to create an inclination-corrected sample of bright exponential galaxies, 
their apparent shapes are confirmed to be consistent with a distribution of mildly 
non-circular disks, with median short-to-long axis ratio 7 ~ 0.22 and median disk 
ellipticity e ^ 0.08. 

Subject headings: galaxies: fundamental parameters — galaxies: photometry — 
galaxies: spiral — galaxies: statistics 
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INTRODUCTION 



Luminous galaxies (Mr < —18 or so) can be coarsely divided into two classes, con- 
ventionally labeled "early-type" and "late-type". Early-type galaxies have redder stellar 
populations and a scarcity of interstellar gas and dust. The majority of luminous early- 
type galaxies are elliptical galaxie s, characterized by srn ooth isophotes and concentrated 
light profiles, well described by a Ide Vaucouleurd (119481 ) profile: log/ oc — r^/^. Highly 



luminous elliptical galaxies tend to be mildly triaxial ellipsoids (as opposed to perfectly 
oblate spheroids); their intrinsic short-to-long axis ratio is typically c/a ~ 0.7 (jRydenlll992l : 
Vincent fc Rvdenlboosh . 



Late-type galaxies have blue stellar populations and relatively large amounts of inter- 
stellar gas and dust. The majority of luminous late-type galaxies are spiral galaxies, char- 
acterized by spiral structure within flattened disks. The disk light profile is generally well 
described by an exponential profile: log J oc — r. Luminous spiral galaxies are mildly elliptical 
(as opposed to perfectly circular) when seen face-on; their intrinsi c short-to-long axis ratio is 



0.25 (Binnev fc de Vaucouleurs 



Lambas et al. 


1992; 


Fasano et al. 


1993; 


Rvden 


2004. 


2006) 



In a color-magnitude (CM) diagram, if the color index is chosen correctly, the early- 
type and late-type galaxies manifest themselves as a "red sequence" and a "blue sequence" , 
respectively. In the Sloan Digita l Sky Survey ( S PSS), the distribution of m — r colors for 
low-redshift galaxies is bimodal; IStrateva et al.l (120011 ) find an optimal color separator of 
u — r = 2.22, when color alone is used as a discriminator between early-type and late-type 
galaxies. Using a f ull CM diagram, the color separator is found to be dependent on Mr 
( iBaldry et al.ll2004l ). with the color separator ranging from u — r 2.3 for galaxies with 
Mr < -21 to M - r ^ 1.8 for galaxies with Mr > -18. 

However, a clean separation between early-type and late-type galaxies using color and 
absolute magnitude information alone is impossible; the red sequence and blue sequence 
overlap in a CM diagram. This overlap results partly from the fact that the color and 
apparent magnitude of spiral galaxies are inclination dependent. Since spiral galaxies are 
disk-shaped and contain dust, an edge-on spiral is both redder and faiii ter than the same 
spiral would be if seen face-on. Thus, as noted by lAlam fc RydenI (120021 ). a sample chosen 
solely by the color criterion u — r > 2.22 will contain dust-reddened edge-on spirals as well 
as intrinsically red ellipticals. 

The overlap between the red sequence and blue sequence would be reduced if we could 
perform an inclination correction on the colors and apparent magnitudes of spiral galaxies; 
that is, if we could convert observed apparent magnitudes into what they would be if the 
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spiral galaxy were face-on. Since pivoting galaxies so that we can see them face-on is an 
impracticable task, we will take a statistical approach to finding the average dimming (AM^) 
and reddening (A(m — r)) of a spiral galaxy as a function of its inclination i. 

In addition to allowing a cleaner separation between early-type and late-type galaxies 
in the CM diagram, a statistical correction for dimming and reddening has other practical 
uses. For instance, a fiux-limited survey will undersample edge-on spirals with respect to 
face-on spirals; a galaxy that is just above the flux limit when it is face on would fall below 
the limit if it were edge on. With a knowledge of AM^ as a function of inclination, it 
is possible to create an inclination-corrected flux-limited sample, by retaining only those 
galaxies that would still be above the flux limit if they were tilted to be seen edge on. 
Since the standard technique for finding the distribution of intrins ic axis ratios of galaxies 



assumes a random distribution of inclinations (see, for instance, (IVincent Sz RydenI l2005l ) 
and references therein), such an inclination-corrected sample is essential for determining 
the true distribution of disk flattening. Without the inclination correction, the scarcity of 
edge-on galaxies would lead to an overestimate of the typical disk thickness. 

In section |21 we describe how we select a sample of SDSS galaxies with exponential 
profiles; this gives us a population of disk-dominated spiral galaxies. The apparent axis 
ratio q of the 25 mag arcsec"^ isophote is chosen as our surrogate for the inclination of a 
galaxy. In section [3l we examine the luminosity function of the SDSS exponential galaxies 
as a function of q. By shifting the luminosity functions until their high-luminosity cutoffs 
align, we can estimate the dimming AM^ of the cutoff as a function of q. By brightening 
each galaxy by the AM^ appropriate to its observed value of q, we can find its approximate 
"face-on" absolute magnitude M/. We then provide linear fits for the mean u — r color as 
a function of M/ for different ranges of q. In section HJ we create an inclination-corrected 
flux-limited sample, which we then use to find the distribution of intrinsic short-to-long 
axis ratios of the SDSS exponential galaxies, confirming that they are, in fact, disks, as our 
analysis assumed all along. Finally, in section [5|, we provide a brief discussion of what the 
form of AMr as a function of q and of A(u — r) as a function of q and M/ implies for the 
properties of spiral galaxies and the dust they contain. 



DATA 



The Sloan D i gital Sky Survey (SDSS) has i maged roughly vr steradians of the sky 



( lYork et al.l I2OOOI : lAdelman-McCarthy et al.l l2008l ). SDSS photome tric data is provided 



in five bands (uq riz) from the near ultraviolet to the near infrared (iFukugita et al.l Il995 



Smith et al.ll2002l ). The SDSS Data Release 6 (DR6) includes 9583 square degrees of photo- 
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metric coverage and 7425 square degrees of spectroscopic coverage (lAdelman-McCarthy et al. 



20081 ). The SDSS photometric data processing pipehne performs a morphological star/galaxy 
separation, with extended objects being labeled "galaxies" and point-like objects being la- 
beled "stars". For each galaxy, in each photometric band, a pair of mod els are fitted to 



the t wo-dimensional galaxy image. One model has a de Vaucouleurs profile (Ide Vaucouleurs 



19481): 



I{R) = h exp (-7.67[(/?//?e)~'/' - 1]) , (1) 

truncated beyond IRe to go smoothly to zero at 8i?e- The other model has an exponential 
profile: 

I{R) = h exp [-1.68(i?/i?e - 1)] , (2) 

truncated beyond 3i?e to go smoothly to zero at AR^. The SDSS DR6 pipeline also takes 
the best-fitting de Vaucouleurs model and exponential model for each galaxy, and finds 
the linear combination of the two that best fits the galaxy image. The fraction of the 
total flux contributed by the de Vaucouleurs component is the parameter fracDeV, which 
is constrained to lie in the interval < fracDeV < 1. The fracDeV parameter thus acts as 
a concentration index, varying from fracDeV = 1 for highly concentrated de Vaucouleurs 
galaxies to fracDeV = for less concentrated exponential galaxies. 

As our measure of the flux of each galaxy, we use the r band "model magnitude"; for 
galaxies with fracDeV > 0.5, the model magnitude is the integrated magnitude of the de 
Vaucouleurs model, and for galaxies with fracDeV < 0.5, it's the integrated magnitude of the 
exponential model. As our measure of the color of each galaxy, we use the u — r color, defined 
as the difference between the u band model magnitude and the r band model magnitude. 
As our measure of the axis ratio of each galaxy, we use the axis ratio of the r band 25 mag 
arcsec"^ isophote. The SDSS DR6 data pipeline finds the best-fitting ellipse to the 25 mag 
arcsec"^ isophote in each band; the semimajor axis and semiminor axis of this isophotal 
ellipse are and -B25, respectively. The isophotal axis ratio ^25 = -825/^25 then provides a 
measure of the apparent galaxy shape at a few times t he effective radius; for g alaxies with 
fracDeV = 0, the average value of A2^/Re is about 2.4 (jvincent &: Rydeii 2005). 



Our full sample of galaxies consists of objects in the SDSS DR6 spectroscopic sample 
that are labeled as galaxies and that have spectroscopic redshifts z > 0.004, to eliminate 
contaminating foreground objects, and z < 0.06, to reduce the possibility of weak lensing 
distortions of apparent shapes, and to eliminate, in practice, the necessity of applying K- 
corrections. To eliminate low-quality redshifts, we require that the SDSS redshift confidence 
parameter have a value zConf > 0.35. To ensure that the galaxies in our full sample are well 
resolved spatially, we require that their photometric data fulfill the criterion r > 6.25rpsf, 
where r is the adaptive second-order moment of the galaxy image and Xpsf is the adaptive 
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second-order moment of the point spread function at the galaxy's location. Our full galaxy 
sample, defined in this way, contains n = 78,230 galaxies. 

Figure [T] shows the color - magnitude diagram for the 78,230 galaxies in our full sample. 
The absolute magnitude Mj. is computed from the model magnitude assuming a Hubble 
constant Hq = 70kms~^Mpc~^ in a flat universe with mass contributing fim,o = 0.3 to the 
density parameter and a cosmological constant contributing = 0.7. No /^-corrections 
are applied to the galaxies in our sample. In Figure [H the division between the blue sequence, 
on the left, and the red sequence, on the right, can be clearly seen. However, the so-called 
"green valley" , between the blue and red sequences, is well populated with galaxies. 

Using the information about surface brightness profiles provided by the fracDeV pa- 
rameter, we can isolate subsamples that are all late-type or all early-type galaxies. Of our 
full sample, n = 36,162 galaxies have fracDeV < 0.1. The CM diagram for this exponen- 
tial subsample is shown in the left panel of Figure [2l The exponential galaxies fall along 
a well-defined blue sequence. Note, though, that the spread in m — r color increases with 
increasing luminosity. For comparison, the right panel of Figure [2] shows the CM diagram 
for the n = 27,567 galaxies in our full sample that have 0.1 < fracDeV < 0.5. These more 
concentrated galaxies fall primarily along the blue sequence; however, there are a significant 
number at the faint end of the the red sequence {Mr ~ —19.5, u — r 2.4). Moving to still 
greater concentration, the left panel of Figure [3] shows the CM diagram for the n = 11,202 
galaxies that have 0.5 < fracDeV < 0.9. These galaxies fall primarily along the red sequence; 
however, there are a significant number at the bright end of the blue sequence {Mr ~ —20, 
u — r ~ 1.6). Finally, the right panel of Figure [3] shows the CM diagram for the relatively 
small number of galaxies in the full sample (n = 3299) that have fracDeV > 0.9. These 
very red, very luminous de Vaucouleurs galaxies represent the high-luminosity end of the red 
sequence. 

Previous studies of the colors, absolute magnitudes, and apparent shapes of late-type 
galaxies in the Sloaii Digital Sky Survey have used different definit i ons of "late-type galaxy" . 



Blanton et al.l (120031 ) . for instance, find that galaxies with ISersid (Il968l ) index n < 1.5 fall 
on the blue sequence of a CM diagram, while those with Sersic index > 3 fall on the red 
sequence. (A galaxy with a perfect Sers i c pro file wi ll have fracDeV 0.3 if = 1.5 and 



fracDeV ~ 0.8 if n = 3.) IChang et al.l (120061 ) and IShao et al.l (120071 ) . in their studies of 



late-type galaxies in the SDSS, find it useful to chose a sample with fracDeV < 0.5. This cut 
in fracDeV, however, does allow some galaxies from the red sequence to enter the sample. 
Since the fracDeV distribution is strongly peaked at fracDeV = 0, we choose to make the 
more stringent cut fracDeV < 0.1 to create our late-type galaxy subsample. The n = 36,162 
galaxies with fracDeV < 0.1, we expect, will provide us with a sample of late- type spirals 
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whose light is dominated by a dusty exponential disk. 



3. 



ANALYSIS 



How do we know that our subsample of (nearly) exponential galaxies are actually flat- 
tened disks? We have, after all, chosen them solely on the basis of their surface brightness 
profiles. Although the disks of spiral galaxies are known to be well fitted by exponential pro- 
files, so are other subclasses of galaxies, such as dwarf ellipticals. We can test the assertion 
that exponential galaxies are fiattened, dust- containing disks by looking at their colors as a 
function of apparent axis ratio. If a galaxy is a disk, then its apparent axis ratio q will be 
smallest when it is edge-on. If the galaxy's disk is dusty, then the galaxy will be most red- 
dened when it is edge-on. Figure H] is a plot of the mean u — r color as a function of absolute 
magnitude Mr for the n = 36,162 galaxies in our exponential subsample (fracDeV < 0.1). 
The galaxies are binned by their apparent axis ratio q. Note that for exponential galaxies 
brighter than Mr ~ —19, there is a noticeable correlation between q and m — r at a given 
absolute magnitude, with smaller values of q corresponding to larger values of m — r. This 
is just what we expect for a population of dusty disk galaxies. However, at Mr > —18, 
there is no correlation between q and u — r. At these low luminosities, the galaxies in the 
exponential subsample are blue (u — r ~ 1.5) dwarf galaxies in which the stars and dust are 
not in orderly thin disks. Thus, in looking for inclination-dependent colors and magnitudes, 
we will look at those high-luminosity {Mr < —19) exponential galaxies for which the color 
actually is inclination-dependent. 

In general, when a spiral galaxy is seen at an arbitrary inclination, it will be fainter at 
visible wavelengths than it would be seen face-on. We may write 



where Mr is the r-band absolute magnitude computed from the actual apparent magnitude 
and redshift, M/ is what the absolute magnitude would be if the galaxy were seen face-on, 
and AMr > is the inclination-related dimming. In general, AMr will be a function of the 
(unknown) inclination i as well as of the detailed properties of the observed galaxy's dust. 
We will assume that the apparent axis ratio q of the 25 mag arcsec"^ isophote will be an 
adequate measure of the inclination. If every spiral galaxy were a perfect oblate spheroid, 
with intrinsic short-to-long axis ratio 7 = c/a, then the inclination i would be uniquely 
determined by the apparent axis ratio q, through the usual relation 



Mr = Ml + AM, 



(3) 




(4) 
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One source of error in this relation is that spiral galaxies don't all have the same intrinsic 
thickness 7. However, an erroneous assumed value for 7 is irrelevant in the case of low 
inclination and high q: a face-on thick disk has the same apparent shape as a face-on thin 
disk. An erroneous assumed value for 7 produces the largest error in i when the apparent 
axis ratio q is small: a circular disk with q = 0.3 is at an inclination cosi = if it's a fat 
disk with 7 = 0.3, but at an inclination cosi = 0.28 if it's an ultrathin disk with 7 = 0.1. 
Another source of error in equation (jl]) is that disks are not perfectly circular. Consider 
an ultrathin disk with e ~ 0.1, a typical ellipticity for a spiral galaxy. When the galaxy is 
face-on, equation (jlj) will yield cosi ~ 0.9, instead of the true value of cosi = 1. When the 
same disk is viewed at a high inclination, the value of cosi we compute will have an error 
~ e cos 20, where (t> is the azi muthal viewing angle measured relative to the intrinsic long 



axis of the disk (iRyded |2006| ). Thus, when we use g as a surrogate for the inclination i, 
we should remember that two spiral galaxies with the same q may have cosz differing by as 
much as ~ e, where e is their average intrinsic disk ellipticity. 

To simplify our analysis, we will start by assuming that AMr is a function only of q, 
and not of M/. That is, we assume that all spiral galaxies with a given q suffer the same 
fractional loss of flux in the r band. This is not necessarily correct - bright spirals may 
well have systematically different dust properties from dimmer spirals, for instance - but it 
provides a place to start. To find AMr(g), we begin by creating a volume-limited subsample 
of our exponential (fracDeV < 0.1) SDSS galaxies. The redshift limit of the volume-limited 
subsample is z = 0.06; the flux limit is assumed to be = 17.77, resulting in a low- 
luminosity cutoff of Mr = —19.4. (The SDSS spectroscopic survey is complete to a limiting 
Petrosian apparent magnitude r = 17.7; for galaxies that are well fitted by the exponential 
model galaxy, the Petrosian apparent magnitudes and the model magnitudes used in our 
analysis have an rms scatter of < 0.1 mag.) We then take the n = 16,363 exponential 
galaxies in the volume-limited subsample and bin them by apparent axis ratio, in bins of 
width Aq = 0.1. The galaxies with 0.9 < q < 1.0 are galaxies that we expect to be nearly 
face-on. 

The cumulative luminosity function for the exponential galaxies with 0.9 < q < 1.0 
is shown as the solid line in the upper panel of Figure O Since we expect AMr = for 
these nearly face-on galaxies, this solid line is our estimator of the normalized luminosity 
function /(M/) for face-on spiral galaxies brighter than M/ = —19.4. For comparison, 
the dashed line in the upper panel of Figure [5] shows the cumulative luminosity function 
for exponential galaxies with 0.2 < q < 0.3; we expect these galaxies to nearly edge-on. 
The luminosity function for the q ~ 0.25 galaxies is different from that of the q ~ 0.95 
galaxies; a Kolmogorov-Smirnov test yields a probability Pks = 4 x 10~^^. This difference 
results from the fact that the luminosity function of spiral galaxies has a cutoff at high 
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luminosity; for edge-on galaxies, this cutoff is dimmed by an amount AM^. Using our 
assumption that dimming (in magnitudes) is independent of M/, we take every galaxy in 
the 0.2 < q < 0.3 and brighten it by the same amount AM^. We then compare this 
artificially brightened luminosity function to the luminosity function for face-on (0.9 < q < 
1.0) galaxies brighter than Mr = —19.4 — AM^,. The comparison is done between the 
cumulative luminosity functions using a Kolmogorov-Smirnov test. The bottom panel of 
Figure [5] shows the comparison between the 0.2 < q < 0.3 galaxies and the 0.9 < q < 
1.0 galaxies, using the optimal shift AM^ = 0.51. The Kolmogorov-Smirnov test yields a 
probability Pks = 0.96 for this shift; our hypothesis of AMr independent of M/ cannot be 
excluded, at least in the luminosity range M/ < —19.91. 

Figure Elgives a summary of the best-fitting absolute magnitude shift AM,, as a function 
of q for the exponential galaxies. Each data point gives the value of AM^ which maximizes 
the Kolmogorov-Smirnov probability when the luminosity function of galaxies in a particular 
range of q are compared to galaxies with q > 0.9. The error bars indicate the range of AMr 
over which Pks > 0.1. Because the high- luminosity cutoff in the luminosity function is not 
an extremely sharp feature, the range of statistically acceptable values of AMr is as large as 
~ 0.4 mag in the low-g bins, where there are relatively few galaxies. A c ommon parameteriza- 



tion o f the AM - q relation, following the Third Reference Catalogue (Ide Vaucouleurs et al. 
199lh . is 

AM = -7Alogg, (5) 

where the value of 7a depends on the waveleng th of observation, and on the population of 
galaxies observed. Ide Vaucouleurs et al.l (Il99ll ) used ■jb = 1-5 for spiral galaxies o f type 
Sc observed in th e B band. For p ure disk systems (type Sc-Sd), iBottinelli et al.l (119951 ) 
found 7b = 1.67. iTully et al.l (119981 ) found dimming that was dependent on luminosity as 
well as on wavelength: for bright spirals in the Ursa Major a nd Pisces clusters , they found 
7b = 1.7, 'jR = 1.3, and 77 = 1.0 in the B, R, and / bands. IShao et al.l (120071 ) . looking at 
fracDeV < 0.5 galaxies in the SDSS, found values of 7 ranging from 7^ = 2.19 in the w-band, 
through 7j. = 1.37 in the r-band, to 7^ = 0.80 in the z-band. 

Our best fit to the logarithmic relation of equation is shown as the dashed line in 
Figure ini the best fitting value of 7 is 7,. = 0.64. Note, however, that the AMr oc \ogq fit is 
not very good. It overestimates the dimming of galaxies with q > 0.4, and underestimates 
the dimming of nearly edge-on galaxies with q < 0.3. In fact, our data are consistent with 
no dimming at all in the r band for galaxies with g > 0.5. A superior fit is provided by the 
solid line in Figure El which represents a dimming proportional to the square of logg: 



AMr = Prilogqf 



(6) 



with Pr = 1.27. Our results are in qualitative agreement with those of iMasters et al.l (120031 ) 
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who find that for galaxies in the Two Micron All-Sky Survey Extended Source Catalog, 
the simple linear extinction law of equation ([5]) gives an inadequate fit to the dimming in 
near-infrared (J, if, and Kg) bands. A better fit is provided by a bilinear function with 
a steeper slope at logg < —0.5 than at log o > 0.5. A good fi t is also provided by a 



quadratic in log q. In the J band, for instance, [Masters et al.l (120031 ) find the best quadratic 
fit is AMj = 0.12 log g -|- 1.14(logg)^. If we attempt a fit of this form to our exponential 
r-band subsample, we find AM^ = 0.25 logg + 1.66(logg)^. However, this does not provide a 
statistically better fit, given the loss of a degree of freedom, than the simpler parabolic form 
of equation Qj. 

Given the dimming correction of equation ([6]), we can compute the "face-on" absolute 
magnitude M/ for every exponential galaxy with measured m^, z, and q: 

M/ = Mr- AMr = Mr- 1.27(logg)2 . (7) 

Figure [7] shows the average u — r color for exponential galaxies, as a function of the corrected 
absolute magnitude M/. We look only at galaxies with M/ < —19, to exclude the dwarf 
galaxies for which the correction is inappropriate. The different colors and line types in 
Figure [7] represent different values of the apparent axis ratio q, just as in Figure HI For 
each range of q, the mean color, (u — r), is linear in the corrected absolute magnitude, 
M/, with brighter galaxies being redder, on average. The q > 0.9 galaxies, which are 
nearly face-on, have a relatively small dependence of average color on absolute magnitude; 
for the q > 0.9 galaxies, increasing the brightness by 1 mag in the r band corresponds to 
reddening the galaxy by ~ 0.1 mag, on average, in u — r. This correlation is a manifestation 
of the dependence of stellar population on galaxy luminosity; in ge neral, more luminous 



spiral galaxies have po pulations that are older and more metal-rich (IBell Sz de Jong 12000 



MacArthur et al.ll2004l ). The nearly edge-on galaxies (g < 0.3), have a greater dependence 
of average color on absolute magnitude; for the q < 0.3 galaxies, an increase in brightness 
by 1 mag in the r band corresponds to a reddening of ~ 0.27 mag, on average, in m — r. 
The steepness of the mean color - luminosity relation for nearly edge-on spiral galaxies 
is a manifestation of the dependence of dust opacity on absolute magnitude; in general, 
more luminous spira l galaxies, being more metal-rich, have greater disk opacity due to dust 



flMasters et al.ll2003[ ) 



By doing a least-squares fit to the function 

{u — r) = a + h{Mf + 20.5) , 



we find the straight lines in Figure [3 the intercepts (a) and slopes (6) are given in Table [TJ 
The mean u — r color for an exponential galaxy with M/ = —20.5 is nearly linear in logg; 
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our best fit is 

a = 1.72 - 0.723 log g . (9) 
The best fitting function of a similar form for 6 as a function of q is 

6 = -0.11 + 0.255 logg . (10) 

We can compute a corrected, "face-on" color, {u — r)-^, for each galaxy using the relation 



iu — rV = (u 



r) + [0.723 - 0.255(M/ + 20.5)] logg , (11) 



where -u — r is the color computed from observations, and M/ is the corrected absolute 
magnitude from equation ([7]). To illustrate the effects of using the "face-on" colors and 
absolute magnitudes on the CM diagram. Figure [8] shows the CM diagram for the exponential 
(fracDeV < 0.1) galaxies in our fiux-limited z < 0.06 sample. The left panel uses the 
corrected absolute magnitudes and colors from equations ([7]) and flTT]) . while the right panel 
uses the uncorrected values. Of the 21,813 galaxies with Mr < —19 plotted in t he right panel 



of Fig ure El 5213 (23.9%) have au — r color redder than the optimal divider of iBaldry et al. 



( I2OO4J ). intended to divide the red sequence from the blue sequence with optimal effectiveness. 



By contrast, of the 23,602 galaxies with M/ < —19 pl otted in the lef t pane l, only 2335 (9.9%) 



have colors redder than the optimal dividing color of iBaldry et al.l (120041 ). 



The color correction of equation ffTTj) is an average correction; the actual correction 
for any individual galaxy will be different, due to the galaxy-to-galaxy variation in dust 
properties. The consequence of the variation is seen in Figure |9l which shows the standard 
deviation in the u — r color, as a function of M/, for different values of q. For nearly face- 
on (g > 0.9) exponential galaxies, the standard deviation in color is luminosity dependent, 

decreasing from cr{u — r) ~ 0.3 at M/ 19 to a{u — r) ~ 0.2 at M/ ~ —21.5. For 

nearly edge-on galaxies (g < 0.3), the standard deviation is greater at any given value of 

M/, ranging from a{u - r) ~ 0.35 at M/ 19 to a{u - r) ~ 0.3 at M/ 21.5. At any 

value of M/, the edge-on galaxies have an excess in a{u — r), caused by variations in dust 
properties, of Act ~ 0.23, added in quadrature. 



4. INTRINSIC SHAPES 

Our analysis has implicitly assumed that bright SDSS galaxies with fracDeV < 0.1 are 
flattened, nearly circular disks. Now that we have an empirical correction for dimming as 
a function of apparent axis ratio (equation [7]) , we can test whether this assumption is self- 
consistent. If a population of galaxies consists of oblate spheroids with a random distribution 
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of inclinations, the observed distribution of appa rent axis ratios, f(q), can be inve rted to 



find the distribution of intrinsic axis ratios, /(7) (lHubblelll926l : ISandage et al.lll970l ) 



For a population of dusty disk galaxies whose apparent magnitude depends on incli- 
nation, selecting a subsample whose inclinations are random requires a little extra care. A 
simple flux-limited subsample will be biased against disks at high inclination, which will 
have a lower flux than a low-inclination disk at the same redshift with the same luminosity. 
Thus, we create a corrected fiux-limited sample; for each SDSS exponential (fracDeV < 0.1) 
galaxy, we ask, not simply whether it is above our flux limit, but whether it would be above 
our flux limit if it were edge-on. If an SDSS exponential galaxy has q = 0.2, we assume that 
it is already edge-on. If it has an observed apparent axis ratio gobs > 0.2, we compute its 
edge-on flux to be 

m^(edge-on) = m,.(observed) + AM^(g = 0.2) - AM,.(gobs) (12) 
= m,(observed) + 1.27[(log 0.2)2 - (log gobs)T ^^3^ 

We create our corrected fiux-limited sample by demanding that this computed edge-on flux 
be greater than = 17.77, the completeness limit of the SDSS spectroscopic survey for 
galaxies. To ensure that our inclination-corrected sample contains only the luminous, nearby 
galaxies for which our AMr{q) correction was computed, we add the addition restrictions 
M/ < -19.4 and z < 0.06. 

The distribution of the apparent axis ratio q for the n = 16,155 exponential galaxies in 
our inclination-corrected sample is shown in the top panel of Figure [TUJ The distribution 



shown as the solid line is e stimated using a nonparametric kernel technique (IVio et al.lll994 



Tremblay &: Merrittlll995l ). To ensure that our estimate of /(g) is smooth, we use a Gaussian 



kernel. T h e ker nel width is chosen using the standard adaptive two-stage estimator of 



Abramsod (119821 ). To ensure that our estimate of /(g) is zero for g < and g > 1, we 
impose reflective boundary conditions at g = and g = 1. The dashed lines in Figure [TD] 
indicate the 98% error intervals on /(g) found by bootstrap resampling. In our bootstrap 
analysis, we randomly selected n = 16,155 values of g, with substitution, from the original 
set of n data points, and then created a new estimate of /(g) for the bootstrapped data. 
After doing 500 such bootstrap estimate, we then determined the 98% error intervals show 
in Figure [TUl At every value of g, 1% of the bootstrap estimates lie above the upper dashed 
line, and 1% lie below the lower dashed line. 

Having a smooth estimate of /(g), the distribution of apparent axis ratios, allows us to 
compute N{'y), the distribution of intrinsic axis ratios, given the assumption that all galaxies 
are oblate or prolate. If the disk galaxies are assumed to be oblate, the relation between 
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/(g) and A^(7) is 



/(?) 



(14) 



where Pohi{(l\l)dq is the conditional probabihty that an oblate spheroid with an intrinsic 
short-to-long axis ratio 7 has an observed apparent axis ratio in the range o — » q -|- dq, 
averaged over all viewing angles. The numerical value of Pobi is ( ISandage et al.lll970l ) 



Pou{q\l) 



(1 _ ^2)1/2(^^2 _ ^2)1/2 



(15) 



if 7 < g < 1, and Pobi = otherwise. Equation (HM is a Volterra equation of the first kind; 
in its discretize d form, it can be i i iverte d by a process of forward substitution to find A^(7) 
given /(g) (see lVincent fc RydenI (120051 ) for numerical details). 



In the bottom panel of Figure [TDl the solid line indicates the estimate of A^(7) found 
by inverting the best estimate of /(g); the dashed lines are the 98% confidence intervals 
found from the inversion of the bootstrap estimates of /(g). The most probable value of 
7, given the oblate hypothesis, is 7 = 0.22. For comparison to these r-band results, an 
inclination-corrected sample of spiral galaxies from the 2MASS Large Galaxy Atla s has a 
most probable thickness of 7^ = 0.17 in the B band and '-fx = 0.25 in the Kg band (iRyden 
20061 ) . A noteworthy property of our estimate of A^(7) is that it is negative for large values 
of 7. The 98% confidence interval is negative for 7 > 0.89; that is, fewer than 1% of our 
bootstrap resamplings give iV(7) > in this interval. This unphysical results permits us to 
reject the hypothesis of perfect oblateness at the 99% (one-sided) confidence level. 

We can approximate a galaxy in our sample not as an oblate spheroid, but as a triaxial 
ellipsoid, with axis lengths a > b > c. The shape can then be expressed in terms of 
two parameters, which we choose to be the short-to-long axis ratio, 7 = c/a, and the disk 
ellipticity, e = 1—6/ a. Once we permit non-zero values of e, we can no longer use the observed 
distribution /(g) to uniquely determine the intrinsic distribution A^(7,e). However, we still 
can do parametric fits to the distribution of shapes. It is found that a useful parameterization 
is A^(7,e) = N^{'j)Ne{e), with A^^ being a Gaussian, 



Ny{'y) (X exp 



(7 - /^t) 
2a^ 



21 



and being a log-normal distribution (lAndersen et al.ll200ll : lRydenll2004l . 120061 ). 



A^. 



e J oc - exp 



(Ine - fif 
2^^ 



(16) 



(17) 



The best values of /i^, a^, /i, and a were determined by a fit to the binned distribution of 
g for the n = 16,155 exponential galaxies in our inclination-corrected sample. The bin width 



- 13 - 



chosen was dq = 0.01. After selecting values of /i^, a^, /i, and 7, we randomly chose n values 
of 7 and of e from the distributions of equation f|T6|) and f|T7|) . For each (7, e) pair, a random 



viewing angle was chosen, and the resulting apparent axis ratio was computed (lBinneylll978l ). 
The model axis ratios were then binned in the same way as the actual, observed axis ratios. 
Repeating this procedure 400 times for each {fj,^,a^, fi,a) set, we calculated the mean and 
standard deviation for the number of model galaxies in each q bin, and computed a 
probability for that particular set of parameters. The best fit we found was /i^ = 0.216 and 
a-y = 0.067 for the distribution of disk thicknesses, with fi = —2.56 and a = 0.91 for the 
distribution of the natural logarithm of the ellipticity. The probability yielded by this set 
of parameters is P ~ 2 x 10~^. In the top panel of Figure [10], the dotted red line shows the 
distribution of q expected from this set of parameters, smoothed with the same kernel width 
as the real data (solid black line). 

The distribution of ellipticity of our best-fitting parametric model, Ine = —2.56 ± 0.91, 
implies a modal ellipticity emod = 0.033, a median ellipticity emed = 0.077 and a mean 
ellipticity eavc = 0.11. This distribution of thicknesses is consistent with the implied ellipticity 
of an inclina tion-corrected sample of spiral galaxies from the 2MASS Large Galaxy Atlas 



(jRyderul2006[ ). if the shape parameter for the 2MASS spirals is the axis ratio of the 25 mag 



arcsec isophote. If the potential eUipticity equals the disk ellipticity, then Ine = —2.56 ± 



0.91 implies roughly one magnitude of scatter in the TuUy-Fisher relation (iTuUy &: Fisher 



19771 ). 



5. DISCUSSION 

We have selected a population of galaxies from the Sloan Digital Sky Survey which 
are at low redshift {z < 0.06), which are relatively luminous [Mr < —19), and which are 
well described by an exponential surface brightness profile (fracDeV < 0.1). As we have 
shown, the properties of these galaxies are consistent with their being a population of slightly 
elliptical disks containing dust. The median dimensionless disk thickness for these galaxies 
in the r band is 7 ~ 0.22; the median disk ellipticity is e ~ 0.08. 

By fitting the luminosity function for galaxies with different apparent axis ratio q, we 
found that the apparent dimming AM^ is not linearly proportional to logg, but instead 
is much better fitted by AMj. oc (logg)^. The dependence of dimming on inclination is a 
valuable clue to the dust properties within disk galaxies. If certain simplifying assumptions 
are made, the expected attenuation as a fun ction of inclination can be computed for model 



galaxies. For instance, iFerrara et al.l (119991 ) assumed that dust had either the extinctio n 



curve found for Milky Way dust or for Small Magellanic Cloud dust (ICordon et al.l 119971 ). 
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They assumed that disks were perfectly axisymmetric, with a horizontal scale length r = 
4kpc that was the same for both stars and dust. The scale height of t he stars was assumed 



to be 2;^ = 0.35, but the scale height of the dust was allowed to vary. iFerrara et al.l (119991 ) 
found that dimming in the B and / bands were proportional to logg only when the dust 
scale hei ght was greater than that of the stars. However, observation of nearby edge-on disk 



galaxies ( Xilouris et al.lE999l ) indicates that the dust scale height is about half the star scale 
height. 



Rocha et al.l (l2008l ). using a Monte-Carlo radiative-transfer code to make calculations of 
internal extinction in dusty galaxies, found that a quadratic dependence of dimming on log q 
provides a good fit for all plausible dust scale heights, scale lengths, and metallicity gradients. 
Since they were using hydrodynamic galaxy models with spiral structure, they were able to 
confirm that nonaxisymmetric structures such as spiral arms did not significantly affect the 
dependence of the total dimming AM on the apparent axis ratio q. 

In our sample of exponential galaxies from the Sloan Digital Sky Survey, once the 
absolute magnitude of a galaxy is corrected for the inclination-dependent dimming, the 
mean u — r color observed is linearly dependent on the corrected M/. For nearly face- 
on galaxies, with q > 0.9, the dependence of m — r on M/ is relatively weak. We find 
b = —0.096; that is, less than 0.1 magnitude of reddening in u — r for each magnitude 
brighter in M/. For the edge-on galaxies, with q < 0.3, the dependence of m — r on M/ 
is much stronger, with b = —0.265. The mean color - absolute magnitude dependence is 
a manifestation of the metallicity - luminosity dependence. High-metallicity galaxies have 
both redder stellar populations and higher dust contents. For face-on exponential galaxies, 
the dust effects are minimized, and we see the effect of metallicity on stellar populations. 
For edge-on exponential galaxies, we see, in addition, the effect of metallicity on the dust 
content. Exponential galaxies with M/ ~ —21.5 have a typical color {u — r) ~ 1.8 when 
seen face-on, placing them at the tip of the blue sequence in a color - magnitude diagram. 
However, the same bright exponential galaxies when seen edge-on will have Mr ~ —21 and 
{u — r) ~ 2.5, a degree of reddening that smuggles them into the red sequence, as usually 
defined. 
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Fig. 1. — Color - magnitude diagram (absolute r band magnitude versus u — r) for all 
galaxies in our flux-limited z < 0.06 sample. The color bar on the right indicates the value 
of logj^Q F corresponding to each color, where F is the fraction of the galaxies in each bin. 
Bin sizes are AM^ = 0.25 and A{u - r) = 0.05. 
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Fig. 3. — Left: As in Figured], but containing only those galaxies with 0.5 < fracDeV < 0.9. 
Right: As in Figure [H but containing only those galaxies with fracDeV > 0.9. 
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Fig. 4. — Average u — r color as a function of for our flux-limited, z < 0.06 sample of 
exponential galaxies. Results are shown for different ranges of q, the apparent axis ratio. 
Error bars represent the estimated error in the mean color. 
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Fig. 5. — Top: Cumulative luminosity function for the galaxies in the volume-limited sample 
with q > 0.9 (solid line) and with 0.2 < q < 0.3 (dotted line). Functions are normalized 
so that F{< Mr) = 1 at Mr = —19.40. Bottom: Comparison of the cumulative luminosity 
functions after the 0.2 < q < 0.3 subsample has been shifted by — AM^ = —0.51. Functions 
are normalized so that F(< Mr) = 1 at = -19.40 - 0.51. 
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Fig. 6. — The points show the best-fitting shift AM^ of the exponential galaxy luminosity 
function as a function of the apparent axis ratio q. The error bars indicate the range of AM^ 
for which Pks > 0.1. The dashed red curve shows AM^ = — 0.64 log g, and the sohd blue 
curve shows AM^ = 1.27(logg)^. 
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Fig. 7. — Average u — r color as a function of the corrected, "face-on" absolute magnitude 
M/ for our flux-limited, z < 0.06 sample of exponential galaxies. Colors and line types are 
the same as in Figure HI Error bars represent the estimated error in the mean color. 
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Fig. 8. — Left: Color - magnitude diagram for fracDeV < 0.1 galaxies, using the corrected 
(m — r)'^ color (equation [H]) and the corrected M/ absolute magnitude (equation [T]) . Right: 
Same as left panel, but using uncorrected u — r color and Mr absolute magnitude. 
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Fig. 9. — Standard deviation in the u — r color as a function of the corrected absolute 
magnitude M/ for our flux-limited, z < 0.06 sample of exponential galaxies. Colors and line 
types are the same as in Figure HI 
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Fig. 10. — Top: Distribution of apparent axis ratios for the inclination- corrected sample 
of SDSS exponential galaxies. The solid line is the best fit: the dashed lines show the 
98% confidence interval estimated from bootstrap resampling. The red dotted line is the 
distribution of q yielded by the best-fitting parametric distribution of elliptical disks (eqns. [16] 
and [T7j) . Bottom: Distribution of intrinsic axis ratios, assuming galaxies are randomly 
inclined oblate spheroids. The solid line is the inversion of the best fit for f{q); the dashed 
lines are the 98% confidence interval from bootstrap resampling. 



Table 1. Mean color vs. face-on luminosity: //(w — r)—a-\- b{M^ + 20.5) 



q 


a 


b 


0.26 


2.15 


-0.265 


0.35 


2.05 


-0.210 


0.45 


1.96 


-0.212 


0.55 


1.90 


-0.184 


0.65 


1.84 


-0.165 


0.75 


1.81 


-0.162 


0.85 


1.77 


-0.127 


0.95 


1.75 


-0.096 



